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We demonstrate theoretically that a charge current and a spin current are generated 
by spin dynamics in the presence of spin-orbit interaction in the perturbative regime. We 
consider a general spin-orbit interaction including the spatially inhomogeneous case. Spin 
current due to spin damping is identified as one origin of generated charge current, but other 
contributions exist, such as the one due to an induced conservative field and the one arising 
from the inhomogeneity of spin-orbit interaction. 
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1. Introduction 

Spin Hall effect 1-6 is one of the most interesting phenomena in spintronics, which enables 
the control of magnetic properties by purely electrical means. The idea is to induce a spin 
current in a transverse direction to an applied electric field by using spin-orbit interaction. As 
an inverse effect, one can expect the conversion of spin current into charge current or electric 
voltage by use of spin-orbit interaction. This effect, namely the inverse spin Hall effect, was 
proposed by Saitoh et al. 7 and indeed observed experimentally in metallic systems and in 
semiconductor (GaAs). 10 One should note, however, that detection of spin Hall effect has so 
far been done by observing magnetization as a result of flow of spin current, and not the spin 
current itself. In the inverse effect, similarly, the electric voltage is measured as a response to 
an time-dependent external field which drives magnetization dynamics. 

The inverse of spin Hall effect was theoretically pointed out by Zhang and Niu 11 and 
Hankiewicz et al., 12 where they discussed a transverse charge current by a gradient of a spin- 
dependent chemical potential (they called this effect the reciprocal spin Hall effect). The 
spin-dependent chemical potential was argued to be related by an optical method. In a junc- 
tion of ferromagnet attached to normal metal, generation of electric voltage by applying an 
alternating magnetic field observed by Costache et al. 13 Theoretical explanation of dc voltage 
was done by Wang et al. 14 as due to the spin accumulation at the interface arising from a 
backflow of pumped spin current into the ferromagnet. Voltage generation from spin dynamics 
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was predicted by Stern in a slightly different context of Faraday's law for a fictitious field of 
Berry's phase. The application to magnetic domain wall was done by Barnes and Maekawa 16 
and Duine. 17 In the Berry's phase mechanism, spin-orbit interaction is not essential but con- 
tributes as correction. 17 

Direct relation between the pumped charge current and the magnetization dynamics was 
investigated recently by Ohe et al. ls They considered a disordered two-dimensional electron 
gas with the Rashba spin-orbit interaction and interacting with dynamical magnetization. 
The charge current was calculated perturbatively. Diffusive electron motion as represented 
by diffusion pole (proportional to 1/q 2 for small momentum transfer, q) was taken account 
of since it leads to logarithmical long-range correlation in two-dimensions. Considering a 
case of uniform Rashba system, they found that pumped charge current had a contribution 
proportional to {S x S), where S is a local spin and (• • •) here denotes average taking account 
the diffusive motion of electrons. The quantity {S x S) represents a spin damping and is 
related phenomenologically to a spin current across the interface in the case of junctions. 19 ' 20 
The result thus supports the idea of inverse spin Hall effect, where charge current is converted 
from the spin current. Qualitatively this contribution to the current was, however, found and 
the inverse spin Hall mechanism turned out to be too naive. Ohe et al. also noted that the 
uniform Rashba system is peculiar, with many cancellations among the Feynman diagrams, 
similarly to the peculiarity known in the spin Hall effect. 21 

In this paper we extend the study by Ohe et a/. 18 to the cases of general spin-orbit 
interaction, including the case of inhomogeneous spin-orbit interaction. The result can thus 
be applied to a finite Rashba system attached to leads. It turns out that such inhomogeneity 
also contributes to a current with different symmetry proportional to (S). We consider a 
three-dimensional case and so do not take account of diffusion ladders which give only small 
contributions unlike in a two-dimensional case considered in ref. 18. 

2. System 

We consider an electron system with spin-orbit interaction and exchange interaction with 
local spins. The local spins can have arbitrary structure and thus we can discuss various 
systems, such as those with one or two ferromagnets attached to a nonmagnet as shown in 
Fig. 1. We consider a disordered system which would be the case of most experiments in 
metallic systems. The total Hamiltonian is H(t) = Hq + H cx (t) + H so + iTi mp , where 
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Fig. 1. (Color online). Two typical systems with spin-orbit interaction and local spins. The arrow is 
a local spin (S) which may have spatial structure and dynamics, and a bottom layer describes 
material with spin-orbit interaction (s-o). 
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Here the annihilation (and creation) operator of conduction electrons in coordinate space is 
ip x (and ipx). The first term describes free electron. The second term H cx denotes the exchange 
interaction, where J ex is a strength of the exchange coupling, S x (t) represents the local spins 
which can have any spatial and slow temporal structure, and er represent Pauli matrices. 
The spin-orbit interaction is represented by H so , where U x is a scalar potential (including a 
factor h 2 1 4m 2 c 2 ) . The last term H imp is the spin- independent impurity scattering which gives 
rise an elastic electron lifetime r = (2irN e n\u 2 / hV) -1 , where u is a strength of the impurity 
scattering, n\ is a number of impurities, iV e is the electron's density of states at Fermi energy, 
and V is system volume. 

3. Charge Current 

The electron velocity operator is defined as v = j^[H, x] ([A, B] represents the commutator 
AB — BA), which reads = — ^^f- + \t\ivry^r<J v ■ The charge current density is defined as 
j c {x,t) = —e^xity^ipxit)), where A^*v*B = (vA^B + A^(vB) and (• • •} is the expectation 
value estimated by the total Hamiltonian H. It is given by 



j ct i(x,t) = eTr 



( d_ d_ 
2m V dx dx' 



dU. 



dx r i 



G < (x,t;x',t) 



(2) 



where Tr{- • •} represents trace over spin indices and G < (x , t; x' , t) represents a lesser Green 
function defined as G < (x, t; x' , t') = ^(V'l'^OV'ccli)) • This charge current satisfies the charge 
continuity equation, 

dp c (x,t 



dt 



+ V- j c (x,t) = 0, 



(3) 



where p c (x,t) (= — e(ifj x (t)ip x (t))) is the charge density. 
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Fig. 2. (Color online). Diagrammatic representations of the charge current at the first order in J cx . 
Dotted lines and wavy lines denote the exchange interaction with local spin (S) and the spin-orbit 
interaction (s-o), respectively. 

Assuming a dirty case (J ox <C h/r and kpU <C H/t, &f being Fermi momentum), we 
carry out a perturbation expansion to calculate the charge current. We treat the exchange 
interaction to the second order and the spin-orbit interaction to the first order. Path ordered 
Green function 22 ' 23 is defined as G(x,t;x',t') = -|{T c {^(t)^(t')}), where T c {- • •} is a 
path ordering operator denned on Keldysh contour C. This Green function satisfies the Dyson 
equation on complex contour, 



where g x (t) denotes free Green function defined as g x - x '(t — t') = — ji(Tc{ip x {t)il> x i(t')}}oi 
where (• • -)o is the expectation value estimated by free Hamiltonian Hq and averaged over 
impurity scatterings. Dyson equation is solved by iteration. If non-interacting, the charge 
current in //-direction is simply proportional to in momentum space. This contribution 
vanishes, for a system is spatial symmetry (we assume this throughout this paper). The 
charge current first order either in the exchange interaction or the spin-orbit interaction also 
vanishes since TV {a} = 0. Therefore the charge current only arises if exchange and spin-orbit 
interactions couple. 
3.1 First order in J ex 

By use of eqs. (2) and (4), contribution from the left diagram in Fig. 2 is given by 




-i V f dt so g x - R (t - t BO ) {[(VrUr) x V H ] • cr} G(R,t so ;x',t'), (4) 



Jc 



ieh 2 J, 



cx 
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where < denotes taking lesser component. Using Tr{<r a <7^} = 28 a p and calculating a lesser 
component, we obtain 

'\x,t) = — — - — e luv > > e 1 "*— -SVfl 



•(Fig.2-left)/ f) _ [eh J cx f 9 _ d \ ^ m dU Rav 

x | [f(uj + n)- f(u)] g^-x^dx-R^+a {J^9R- x >,u>+^j 

- f(uj + tyg x _ Xu ,g X -R,w+n [-Q^9 T R- x i,u>+n 

where «J*) = ^»^Y M (<7fc,u, = G7fc,J* = - £ fc + £ F + f where 

£fc = h 2 k 2 /2m and £f is Fermi energy), /(w) is the Fermi distribution function which is given 
as f{oj) = 0(— uj) at zero temperature (0(uS) is a step function), and Sx,n is Fourier transform 
of Sx(t)- Assuming that local spins vary slowly (Q <C t -1 ), we obtain 
.(Fi g .2-icftw f) „ eh 2 J cx fd^ 8\ f d \ 

\ / /a; —x 

(7) 

where ^ = g X)U =o- We neglected terms containing only g T, s or <7 a 's since they are higher order 
of h/e-pT < 1 as compared with mixed terms. Other contributions are similarly calculated 
and the whole contribution in Fig. 2 is obtained as 

j£^\x,t) = ^e,^ReJ2Sx(t)9 x -x9x- x 



+ 



eh 2 J cx d 



-e 7 ^Re 2^ -Q^^xK^gx-x \-q^9x-r J gR- x 

.X", -R 

w<« £ (sir &-«) (4:*-) ■ (8) 



X,R 



The first term corresponds to the right diagram in Fig. 2, i.e., correction of the current vertex 
due to the spin-orbit interaction. In momentum space, eq. (8) reads 

xRe£ [(fe-|)^F^_ q5 ^^ +p -^)+^^_ q « + J , (9) 
where S,(t) = £ Ex^' X ^(t) and tf p = f Zr^ R Ur. 
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Fig. 3. (Color online). Charge current at the second order in J cx . (a) and (b) are in the presence and 
the absence of spin-orbit interaction, respectively. 



3.2 Second order in J cx 

Contribution from the second order in J ex shown in Fig. 3(a) is calculated similarly by 
using Tr{er Q <7^cr 7 } = i2e Q/ 3 7 as 



2eJ, 



ex ^ x 



TT 



Im [s Xl (t) x S X2 (t) 



9x-x L 9xi - x 2 9x 2 -x 



Xi,X2 



eh 2 J 2 d v-^ dU R 

+ — Im E 



irm 



dx u ^ dR r > 

M X U X 2 ,R 



S Xl (t) x S X2 (t) 



2-X 



Im 



dU 



R 



X lt X 2 ,R 



S Xx {t) x Sx 2 (t) 



+9 X -R (m^Ir-Xi) 9x 1 -x 2 9x 

+9 x -x 2 9x 2 -R idTv^R-Xx) 9x x -x 

+ ( k ~£~9 x -r) (dW9 T R-x 1 ) 9x 1 ~x 2 9x 2 -x 
+ ( K m~9x-x 1 ^j 9x 1 -x 2 9x 2 -R {dW9R-x) 
, + \£^,9 T x -x 2 ) 9x 2 -R (m^r-xJ 9^-x 



In momentum space, this is written as 
\2eh 2 Jl 

Q,Q,P 

Q+p 



jS is - 3a) 0M) 



irmV 



eyxr, E e- i(Q+p) 'W P \s Q - q (t) x S q (t) 



Im E 



(ii) 



+ [k + ^^{k + qf glg% +q [g% +q+p g% +Q+p + g% +q g% +Q 
_ ( fc + 2+2^ (fc + Q)^ +g ^ +Q (^ +Q+p - ^ +Q ) 
The charge current arises even without the spin-orbit interaction at the second order in 
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J cx . This contribution is shown in Fig. 3(b) and given as 



X\,X 2 

2eh 2 Jl 



9x-x 1 9x 1 -x 2 9x 2 -x 



„ n Im E [^W.^Wl^-x.&.-x^^-,). (12) 

This current is simply due to chemical potential shift by the exchange interaction. It contains 
a scalar product, S ■ S, in contrast with the contribution of Fig. 3(a), and so is small if local 
spins are spatially slowly varying. In momentum space, eq. (12) is written as 
i2eh 2 Jl 



i(Fig.3b), 

JCfl 



( *' t] = £ e ~ lQX i 8 *-*® ■ S q (t)] Im Wfc + f) 9l9% +q 9% +Q . (13) 

q,Q k ^ 



The total charge current is given as j c = j^ lg ' 2 ^ + j'^ Flg ' 3a ^ + ^( Fl s- 3b ) _ xhis current contains 
a component which originates from a conservative scalar potential (or electric field) created 
by the spin dynamics and the spin-orbit interaction, such as the second terms of eqs. (8) and 
(10) and the first term of eq. (12). This scalar potential (or electric field) is fictitious, which 
only acts if the charge current has spin degrees of freedom. The existence of this conservative 
part is in contrast to the adiabatic case in ref. 16. 

The above results are quite similar to those in two-dimensions. 18 The difference is that 
averaging over long-range diffusion motion in two-dimensions is replaced in three-dimensions 
by a short-ranged average within the scale of electron mean free path. 

4. Spin Current 

We consider the spin current arising from exchange interaction. We neglect here spin-orbit 
interaction because we are interested in conversion of spin current into charge current by spin- 
orbit interaction. Spin current density is defined as j^(x,t) = f(V4(£) ^ir ^ i'xif)) {{A-,B} 
represents the anticommutator AB + BA and = — ^gf-), ^ e -' 

h 3 ( d d 



(14) 



In contrast to the charge continuity equation (eq. (3)), this spin current satisfies the equation, 



dt 



+ V-3 v B {x,t) = 1?{x,t), (15) 



where p v s {x, t) (= f is the spin density and T s v [x, t) (= ^Ut)[H ex ,a"]ip x {t))) 

is the torque density (or spin source) due to exchange interaction. Spin is thus not conserved. 
(Definition of spin current has ambiguity when spin-orbit interaction is taken account, 24 ' 25 
but the continuity equation, eq. (15), is always satisfied with T s redefined properly.) 

Diagrammatic representations of the spin current at the first and the second order in J cx 
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Fig. 4. (Color online). Diagrams of the spin current, (a) and (b) denote contributions from the first 
and the second order in J ox , respectively. 



are shown in Fig. 4. These contributions are given by 

h 3 J, 



4 Fig - 4a) (^) 



2imi 



ex ImV^(% 



x 



x-X 



_d_ 

dx. 



(16) 



•KFig.4b) , 
Js/J ' 



X,t) 



lixm dx, 



Re £ [s Xl (t) x S X2 (t) 



9x-x 1 9x 1 -x 2 9x 2 -x 



X\,X 2 



irm 



X\,X2 

They are written in momentum space as 



Re [Sxiit) x Sx a (t) 



9x-x 1 9x 1 -x 2 



dx, 



-9x, 



2-X 



■i/(Fig.4a) 

JSfM 



(x,t) 

ih 3 J 2 



2-nmV 



i : \ (; \Q :r. 



irmV ^— ' 



S Q - q (t) x S q (t)YReY, (k + 



Q 



(17) 



(18) 



9k9k+ q 9k+Q- 



(19) 

The total spin current contains terms like (Sxi) an d (&Xi x &x 2 ), where (• • •) denotes average 
over position X\ and X2 within the length scale of electron mean free path. Comparing these 
expressions with those of charge currents (eqs. (8) and (10)), we see that they are related, 
suggesting the inverse spin Hall effect. In fact, the first term in eq. (8) (eq. (10)) arises from 
(S) ((S x S)), similarly to the spin current of eq. (16) (eq. (17)). These charge currents have 
extra factor of VLXc which justifies the idea of inverse spin Hall effect, conversion of spin 
current into charge current by spin-orbit field, E so = — VLXj,. However, other (second and 
third) terms of eqs. (8) and (10) do not have direct correspondence to spin current. 



5. Smooth Spin Structures 

Let us first consider a case U p is spatially smooth, i.e., p <C £ , where I is the electron 
mean free path. We obtain the charge currents as 

q k 

(20) 
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: ^'"^Y.^ ' Q " [SQ- q (t) x S q (t) 

Q,Q 



nrnV 



<J,Q 



s 2 ^ 



W (k + q) A + (k + ^ fc A ^) 2 (ff| +9 ) 2 ^+Q 

fc + f )^ [k^Q x - k^] (g k ) 2 g% +q (g% +Q ) 2 

+ (fc + [ (fc + g)7 (fc + Q)A + (fe + qV {k + Q) Aj ^ (5 a +g)2(5 a +Q) 2 

+ 0(p 3 ). (21) 

We also assume local spins vary slowly in space, and explore behavior of the currents in detail. 
Carrying out q and Q expansion assuming as q, Q <C £ _1 , the currents are given as 



iefi 2 J< 

3irmV ^ Pri dx^dx x 



d 2 IJ 

cx u u a 



q k 



k 



3irV ^" dxv 



q,Q 
l2eh 2 J 2 



ImA 



M7A 
Q,Q 



TimV 



Q,Q 



J2e~ iQx S Q - q (t) • S q (t) ImB^ Q , (22) 



;fc3 7 



SirmV 



^J 2 X 
irmV 



<J,Q 



where 



4M7A 



<W (9 + Q) A + <W (9 + <5) 7 J E e ^(<?|) 4 

fc 

<W (Q + Q) A + S, x (q - QV] "£e k (g k ) 2 (g%f 

k 

(q + Q) X + 5, x (q + Q^l E e 2 k gl(g%f 



+ 



's m (2q + Q) x + ^ (2g + QV] E 4(£) W> 
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2 3 

Here we used {k a k@) = 8 a pk 2 /3 where the average is over direction. Sums over k are carried 
out as 

fc 



^3 



^ 3 i3vriV c r 2 27r7V c e F r 3 

fc 

i2vriV ceFT 4 



h 3 



fc 
fc 

fc 



2^ 4 ' 

i57riV c eFr 4 
i2vriV c r 2 



h 2 ' 
37riV c r i2vriV c e F r 2 



2fr 



h 2 



where N e is the density of states given as mVk F /2-K h . We finally obtain 



jcfj,( x i t) 



AeN c J cx e F T 3 
3hmV 



d 2 U x dS x (t) d 2 U x dS x (t) 



^8x^8x1 dxi e ^ v dxvdx* dx x 



AeN c J 2 x r 2 3U X 



S x (t) x S x (t) 



8eN c J 2 x e F T 4 
9h 2 mV 



+ 9 

v 

+ ! 



2eiV e J c 2 x r 2 

my 



<9x. 



(24) 



2hN c J cx e F T 2 dS£(t) 
3mV dx„ 



4N e J 2 x e F T 3 [ 3 



3mV I <9x 



X S^t) 



+ 



0x,< 



(25) 



We see in eq. (24) that charge current driven by uniform spin-orbit field (_E so = — W x ) is 
perpendicular to spin damping, S x S, while discontinuity of E so induces new S components. 
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In a case that local spins are spatially uniform, S x (t) = S(t), the spin current totally 
vanishes and the charge current only arises as 



\x,t) = - 



ieN.Jlr 2 dU x 



S(t) x S(t) 



(26) 



5.1 Rashba system 

We apply above result to the case of Rashba spin-orbit interaction. Rashba spin-orbit 
interaction can be realized in various three-dimensional systems without inversion symme- 
try. 26,27 Assuming uniform spin-orbit field in z-direction, —'VU X = E so e z , we obtain the 
charge current as 

OS x {l.) 



J 



Rashba ^ ^ _ jeN c E so J 2 X T 2 



c/i 



h 2 v 



fll/Z 



S x (t) x S x (t) 



" _ 2eN c J 2 x r 2 
mV 



S x (t) ■ 



dx, 



(27) 



It is interesting to note that in the case of inhomogeneous spin-orbit field, an additional current 
arises as 

4e7V e J C x£FT 3 



3hmV 



-flUZ 



dE so ds»(t ) 

dx x dx x 



— e 



■yuz 



dE so ds&t) 

dx^ dx"i 



+ 



8eN c J 2 x e F r 
9h 2 mV 



, 8E so J a d 

"l-e/ii/2 dx x \^ 9x \ 

, dE so f -I q d 



S x (t) x S x (t) 
S x {t) x S x (t) 



+ 9 

V 

+ 9 



S x {t) x ~~ 



8x1 



'} 



(28) 



This additional current appears, for instance, at the interface between Rashba region and 
leads when normal leads are attached. 

6. Equilibrium Components 

In the above discussion on pumped currents, terms like (g T ) n and (g a ) n were negligibly 
small. These contributions in fact correspond to the equilibrium contributions which exist 
even if local spins have no dynamics, S x (t) = S x . The currents are then given by 

f)TT 

j go {X ) = 4e j exe ^^im £ s x £ /!-•;/,;; x .jfx ,., 



Aeh 2 J ex 

H e 

m 



_d 

dx 



8eh 2 J e 



m 



Im E ~Q]^ S x'^2f( u; )g x -x,w (-q^9x-r,uj) 9r- x ,u, 

X,R to ^ ' 

9Ur sr^ t( x a f ® a \ f 3 a 



(29) 



4 cq) w = ^7^- Im E [^^^^/(^/w^ 



M Xi,X 2 



m 



Im J] x Sx 2 ]"X)/H^-Xi, w 55ci-x 2) a, (^4-,,J • 



(30) 
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The equilibrium charge current thus does not contain (Sxi x Sx 2 ) term (ji° q ^ oc (Sx)), while 
spin current does not have (Sx) term (ji eq ^ oc (Sx-i x Sx 2 ))- I n momentum space they are 
written as 

$°(*) = -^^^E e " i(9+P) 'V^5:/(a;) 

q,p <x> 



xIm X] _ f ) ^' g k~q,Lo9k,w(9k+p,w - 9k,w) + g%-q iUj g%,Lo9k+p,u 

k ^ 



(31) 



q,Q w k 'A 4 

(32) 

The equilibrium currents calculated above correspond to an initial response to local spins 
suddenly attached to the system. In reality, as soon as the equilibrium charge and spin currents 
appear, the spin configuration starts to evolve and dynamical currents appear. 

In a spatially smooth local spins and spin-orbit interaction, the equilibrium components 
of the charge current contain at least three derivatives, (VC/a;)(V 2 S x ) or (V 2 [/ a; )(VS' cc ) or 
(V 3 U X )S X , and are negligibly small. The equilibrium spin current, in contrast, is large and 
given as 



h 2 N I 2 

Js ^ [X) 12mVe F 



S x x 



dS x 



1 v 



(33) 



This result indicates that the equilibrium component of spin current cannot be converted into 
charge one. This is physically reasonable since equilibrium spin current can be dissipationless 
while charge current always has dissipation. Dynamical spins are thus essential in pumping of 
charge current. 

7. Conclusion 

We have discussed theoretically that charge and spin currents are induced by spin dynam- 
ics in the presence of spin-orbit interaction. The spin current pumped from dynamical spins 
was found to contain two contributions, (S) and (S x S), consistent with phenomenological 
argument. 7 ' 19 We found that pumped charge current has a component flowing perpendicu- 
lar to both spin-orbit field (E so ) and average spin polarization (S x S). This contribution 
is thus explained as due to conversion of spin current (oc (S x S)) into charge current, i.e., 
the inverse spin Hall effect. For a pumping mechanism by use of local spins, dynamical spins 
are essential because the static spin current is not converted into charge current. The charge 
current has another component originating from a fictitious conservative field which only acts 
if the charge current has spin degrees of freedom. This current is not associated directly with 
spin current. We also found a novel charge current proportional to dE so , arising from inho- 
mogeneity of spin-orbit field. This contribution would be essential in actual experiments on 
finite size spin-orbit system attached to leads. 
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